Negative refraction in (bi)-isotropic periodic arrangements of chiral SRRs 
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Bi-isotropic and isotropic negative refractive index (NRI) 3D metamaterials made from periodic 
arrangements of chiral split ring resonators (SRRs) are proposed and demonstrated. An analytical 
theory for the characterization and design of these metamaterials is provided and validated by 
careful full-wave electromagnetic simulations. The reported results are expected to pave the way to 
the design of practical 3D bi-isotropic and isotropic NRI metamaterials made from a single kind of 
inclusions. 
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Most proposals for developing negative refractive index 
(NRI) metamaterials (e.g. qJ, Q, make use of two 
kind of inclusions, one of them to provide negative per- 
mittivity and the other to provide negative permeability. 
On the other hand, artificial media made from metallic 
chiral inclusions have been known since far by physicists 
[3| , Q. It is known that they present the most general 
linear and reciprocal constitutive relations [y| which, as- 
suming 3D spatial isotropy and a time dependence of the 
kind exp(— juit), can be written as 

(1) 
(2) 



D = E (l + Xe)E + jV £ oMo kH 

B = -j^/soVo «E + /Lt (l + Xm)H , 



where Xe, Xm and K are the electric, magnetic and cross- 
polarizabilities of the metamaterial which, for lossless 
media, are real quantities. Media satisfying (fTJ) and ([2]) 
are usually called bi-isotropc, to differentiate from usual 
isotropic media with n = 0. Since they present simultane- 
ous electric and magnetic properties, it has been guessed 
that they may open the way to NRI metamaterials made 
of a single kind of inclusions. To the best of our knowl- 
edge, the first proposal in this direction was made by 
S.Tretyakov in and further developed in [§]. Also, in 
[t| a mixture of metallic chiral inclusions and wires was 
proposed in order to obtain negative refraction. 

Eigenwaves propagating through reciprocal bi-isotro- 
pic media are right- and left-handed circulary polarized 
plane waves, whose dispersion relation is given by [6| 



k = ko (y/fi r e r ± k) 



(3) 



(1 



where ko — w^/eoMo, and e r — (1 + Xe), 
Xm)- In order to reduce forbidden bands of propagation 
coming from complex values of k^, it is desirable that 
X e (tj) ~ Xm(u) so that [i r and e r always have the same 
sign. As explained in Q and [nj, usually this condition 
also implies that 



(4) 



The conditions for negative refraction in tr ansp arent bi- 
isotropic chiral media were investigated in [111 ]. In loss- 



less bi-isotropic media satisfying the balance condition 
(|H) , these conditions reduces to e r < 0.5 and /j, r < 0.5 
[lfj, which is less restrictive than the usual condition, 
e r < and \x r < 0, for conventional isotropic medium 
(i.e., a medium with k, = 0). The price to pay for this 
enhanced bandwidth is that only one of the eigenwaves 
of ([3]) exhibits negative refraction, the other presenting 
positive refraction [loj . 

Most previously proposed and studied bi-isotropic ar- 
tificial media (see, for instance, 0], [lj| and references 
therein) were random arrangements of chiral inclusions. 
Random arrangements have the advantage of ensuring 
isotropy provided the number of inclusions per unit vol- 
ume is high enough. However, multiple scattering in ran- 
dom arrangements of inclusions affects the effective con- 
stitutive parameters of the mixture, leading to an imag- 
inary part of the susceptibilities to account for radiation 
losses @. This important drawback of radiation losses 
can be avoided by using periodic arrangements of chi- 
ral scatterers. Periodic arrangements also provide higher 
density of scatterers, which will result in stronger electro- 
magnetic responses and wider NRI frequency bands. In 
addition, periodic arrangements can provide much better 
realiability and reproducibility than random mixtures. 
Finally, wave propagation through periodic structures 
can be studied by means of standard electromagnetic 
solvers, thus providing a validation method for the ho- 
mogenization procedure that has no counterpart in ran- 
dom structures. For all these reasons periodic arrange- 
ments of chiral inclusions are preferred over random mix- 
tures. However, care must be taken to obtain a truly 
isotropic periodic metamaterial in order to keep the ap- 
propriate symmetries of the crystal; otherwise couplings 
between elements may destroy isotropy even if the struc- 
ture appears as isotropic in the dipolar approximation 
In it has been shown that a cubic arrangement 
of SRRs satisfying the thetraedron or (in Schoenflies no- 
tation) T-group of symmetry provides an useful basic de- 
sign for this purpose. Fig.l.a shows a cubic arrangement 
of chiral SRRs satisfying this symmetry. It is a combi- 
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FIG. 1: (a) Cubic arrangement of chiral SRRs satisfying T- 
group of symmetry (namely, any of the inbody diagonals is 
a third order rotation symmetry axis), and (b) detail of the 
chiral SRR design. 



nation of six identical chiral SRRs, which are shown in 
detail in Fig.l.b. The proposed chiral SRR is a modi- 
fication of the original SRR design [l3j that substitutes 
the edge-coupling between rings by a broad-side coupling 
14j . also introducing some helicity. It should be noticed 
that the SRRs on opposite sides of the cube of Fig.l.a 
are identical, which makes it possible to develop a peri- 
odic metamaterial with a single cubic lattice from this 
design. The unit cell of this metamaterial has periodic- 
ity a (the edge of the cube) and will contain three chiral 
SRRs. Note that the proposed chiral-SRRs is expected 
to be a feasible "practical" design because of its possi- 
ble manufacture via standard planar-circuit fabrication 
techniques 

The chiral SRR polarizabilities are defined by 
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(5) 
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where (-) T means transpose, p and m are the electric and 
magnetic dipoles induced in the chiral SRR, Ej and Bi 
are the electric and magnetic local fields seen by the SRR, 
and Onsager symmetries (IBj for the magneto-electric 
susceptibilities a = — (a ) T have been explicitly in- 
troduced. Analytical expressions for these polarizabili- 
ties can be obtaine d by a straighforward extension of the 
theory reported in Q . The results are [l6| 
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(7) 
(8) 
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where 



is the external radius of the inclusion, r = 
r ex t — w/2 is the average radius, and X(lo) is a common 
resonant factor given by 



u 2 



lu 2 + jujR/L 



(11) 




FIG. 2: Theoretical dispersion curves for the eigenwaves of 
the periodic metamaterial made of the periodic repetition of 
the cubic arrangement of Fig.l.a, forming a single cubic (sc) 
lattice with peridocity a. SRR parameters are a/r cxt = 2.87, 
w/r cxt ~ 0.5, and t/r cyi t — 0.47. Dispersion curves are drawn 
inside the first Brillouin zone of the structure. 



The frequency of resonance ujq = 1/v ' LC of the inclu- 
sion, as well as the self inductance, L, the capacitance, 
C, and the resitance, R, coincide with those of the afore- 

From 0- 



14] 



mentioned broadside-coupled SRR 
the average magnetic-, electric- and cross-polarizability 
per particle in the isotropic metamaterial are easily ob- 



a, 



73, 



a: 



= /3 + a-/3 + cC/3, 



tained as {a n 

and (a cm ) = cv™/3 respectively. From these expre- 
sions, the susceptibilities of the metamaterial can be ob- 
tained from a straighforward extension of Lorentz local 
field theory that takes into account the presence of cross- 
polarizabilities. This theory leads to the following equa- 
tions for the macroscopic polarization, P, and magneti- 
zation, M, of the metamaterial: 



P = A' i K • y^-) • /io(a cm ) + y 



M = -V { /,„,„„,} ( H + y^j - (a cm ) + ^ 



(12) 

} 

(13) 



where N = 3/a 3 is the number of particles per unit 
volume, and E, H are the macroscopic fields. Us- 
ing (fT2|) -([T3 | the metamaterial susceptibilities are di- 
rectly obtained from their definitions %e = E/(£nP) and 
Xm = H/M. The balance condition ^ is translated 
to the chiral-SRR polarizabilities as (a m ) = |(a cm )|/c = 
(a c ) I c 2 where c is the velocity of light. If the effect of the 



non-resonant polarizability (|T0|) is neglected, the above 
condition reduces to the following approximate balance 
condition: 

t\ = (rrr) 2 , (14) 
where Aq is the wavelength at resonance. 
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FIG. 3: Band structure and dispersion relation for the meta- 
material of Fig. 2 obtained from full-wave simulations along 
the path R-r-X-M-R. 
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FIG. 4: Dispersion curves along the F-X, T-M, and F-R di- 
rections for the eigenwaves of the balanced metamaterial of 
Figs. 2 and 3. 



In the following a particular example is analyzed in or- 
der to illustrate and validate our analytical theory. The 
chosen metamaterial parameters satisfy Eq. [TJ] and are 
given in the caption of Fig. 2. This figure shows the the- 
oretical dispersion curves for the eigenwaves of the bi- 
isotropic metamaterial. A region of backward- wave prop- 
agation for one of the eigenwaves can be clearly observed. 
This frequency band corresponds to the aforementioned 
condition e r , [i r < 0.5. Inside this band there is also a 
small forbidden band gap (see the gray region in the inset 
of the figure) that corresponds to complex values of the 
propagation constants (J3j) . This small band gap is due 
to the aforementioned approximations implicit in Eq. 1141 
The straight horizontal line at the frequency where the 
propagation constant of one of the eigenwaves vanishes 
(corresponding to the condition e r [i r = k 2 ) represents a 
fully degenerate longitudinal wave with k x E = and 
kxH=fl. 

The same structure has been simulated by using CST 
Microwave Studio and its band structure is shown in 
Fig. 3. A good qualitative agreement is found between 
Figs. 2 and 3. In both figures a frequency band of 
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FIG. 5: Dispersion curves along the F-X direction for the 
transverse eigenwaves of two unbalanced metamaterials. Pa- 
rameters are tha same as in Figs. 2 to 4, except for t/r cxt , 
which is taken as t/r cxt = 0.3 (dashed lines) and t/r cxt = 0.1 
(solid lines) 
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FIG. 6: Theoretical dispersion relation for plane waves prop- 
agating in the racemic metamaterial with the single cubic 
lattice shown in the inset. Inset: Racemic periodic structure 
made of two interleaved single cubic lattices of cubes as that 
of Fig.l.a. White cubes are identical to that of Fig.l.a, and 
gray cubes are similar cubes made of chiral SRRs of oppo- 
site handedness. Chiral SRR parameters are as in Fig. 2 and 
periodicity of the structure is 2a. The fully degenerate longi- 
tudinal mode is now located at e r ^r = 

backward-wave propagation is observed for one of the 
eigenwaves. Also a small frequency stopband appears 
in Fig. 3 inside this backward-wave passband (because of 
the approximate balance conditition employed) in agree- 
ment with the theoretical predictions of Fig. 2. In order to 
show the isotropy of the structure, the dispersion curves 
along the T-X, T-M and T-R contours are shown in Fig. 4. 
The isotropy of the dispersion relation for the transverse 
modes becomes apparent from these curves, except for 
high values of the propagation constant, where spatial 
dispersion should affect the dispersion relation. The lon- 
gitudinal mode of Fig. 2 also appears in Figs. 3 and 4, 
although now spatial dispersion destroys degeneracy and 
isotropy. It is illustrative to compare the band struc- 
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ture of the reported balanced structure with those of an 
unbalanced one. Thus, Fig. 5 shows the dispersion dia- 
grams for the transverse modes of two unbalanced struc- 
tres (which were obtained by varying t in Fig.l). It can 
be observed how, as t decreases, the forbidden frequency 
band increases until it occupies all the backward-wave 
region, giving a behavior similar to that of a negative- 
/j, split ring metamaterial [13] • This is expected, since 
the electric and magneto-electric polarizabilities of the 
chiral-SRR ©-© decrease with t, until the chiral-SRR 
becomes an almost purely magnetic inclusion. Fig. 5 illus- 
trates the tolerance of the proposed design to deviations 
from the balance condition (fl4|) . As it can be seen in the 
figure (dashed lines) this tolerance is quite good. 

Let us now explore the possibility of designing a 
racemic mixture of chiral SRRs while keeping the neces- 
sary symmetries to ensure an isotropic behavior for the 
metamaterial. For this purpose let us consider a single 
cubic lattice of peridicity 2a made of cubes as that shown 
in Fig.l. a, and another similar cubic lattice made of the 
same cubes but of opposite handedness Both cubic lat- 
tices can be interleaved as shown in the inset of Fig. 6 
to give a single cubic lattice whose unit cell is formed by 
two cubes of opposite handedness. Since the whole struc- 
ture is symmetric after inversion (the centers of symme- 
try are the corners of the cubes), the cross suceptibility 
k (which is a pseudo-scalar) must vanish [l^|. Therefore 
the whole structure is a balanced isotropic metamate- 
rial with k = 0. Fig. 6 shows the theoretical dispersion 
relation (k = fco^/e r /i r ) of the proposed structure. A fre- 
quency band of backward- wave propagation, correspond- 
ing to the condition e r , fx r < 0, can be observed in the 
figure. 

When designing a NRI metamaterial it is not sufficient 
to obtain some amount of negative refraction but it is also 
required to have small reflectance so that the refracted 
beam has a significant amplitude. As it is well known, in 
isotropic media the reflectance is governed by the wave 
impedance Z = Zq fi r /e r , where Zq = \/JioJeo * 377 f2 
is the wave impedance of free space. For bi-isotropic me- 
dia the same expression Z = Zq^J Jirjsr holds [6]. Since 
for balanced metamaterials it is imposed that Xe — Xm 
in the frequency range of interest, it directly follows that 
Z(uj) ~ Zq(u>), which ensures good matching to free 
space. This last feature is an additional relevant advan- 
tage of the proposed balanced bi-isotropic and isotropic 
NRI metamaterials. 

In summary, periodic NRI metamaterials based on chi- 
ral SRRs have been proposed and demonstrated. An 
analytical theory has been developed for the character- 
ization of such metamaterias, which has been validated 
by careful full wave electromagnetic simulations. It has 
been shown that single cubic lattices of chiral-SRRs can 
provide bi-isotropic NRI metamaterials with a well de- 



fined frequency band of backward-wave propagation for 
one of its plane-wave eigenstates. Besides, a racemic sin- 
gle cubic lattice of chiral SRRs satisfying the apropriate 
symmetry (T^ group in Schoenflies notation) has been 
proposed in order to provide NRI isotropic metamateri- 
als whith a well defined band of backward wave propaga- 
tion for all its plane-wave eigenstates. In both cases good 
matching to free space is expected at all frequencies. 
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